Rotating inclined cylinder and the effect of the tilt angle on vortices 
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(Dated: June 4, 2009) 

We study numerically some possible vortex configurations in a rotating cylinder that is tilted with 
respect to the rotation axis and where difi'erent numbers of vortices can be present at given rotation 
velocity. In a long cylinder at small tilt angles the vortices tend to align along the cylinder axis 
and not along the rotation axis. We also show that the axial flow along the cylinder axis, caused 
by the tilt, will result in the Ostermeier-Glaberson instability above some critical tilt angle. When 
the vortices become unstable the final state often appears to be a dynamical steady state, which 
may contain turbulent regions where new vortices are constantly created. These new vortices push 
other vortices in regions with laminar flow towards the top and bottom ends of the cylinder where 
they finally annihilate. Experimentally the inclined cylinder could be a convenient environment to 
create long lasting turbulence with a polarization which can be adjusted with the tilt angle. 

PACS numbers: 67.30.he,47.32.-y,67.25.dk 
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I. INTRODUCTION 



Under uniform rotation the equilibrium state of a su- 
perfluid consists of a regular array of rectilinear quan- 
tized vortices along the rotation axis. This is strictly 
true only when the top and bottom plates of the rotat- 
ing container are perpendicular to the rotation axis and 
its cross section has fixed translationally invariant form 
along the rotation axis. Near the walls that are parallel 
to the rotation axis the vortex array is surrounded by a 
small vortex free layer of thickness that is of the order of 
the vortex spacing 1, 2]. 

The usual setup to study vortices in rotation is to use 
a cylindrical container and to assume that the cylinder 
axis is aligned exactly parallel to the rotation axis. Nev- 
ertheless, experimentally the rotation axis is always tilted 
by some small unknown angle owing to uncontrolled me- 
chanical misalignment. Even if the misalignment is small, 
one is tempted to ask the question: What is the effect 
of the tilt and how does it affect the steady state vortex 
configuration? Fortunately, it turns out that cylindrical 
symmetry is not immediately broken at any arbitrarily 
small tilt angle. A sensitive test of the tilt angle is the 
measurement of the number and configuration of vortices 
in the equilibrium state of rotation [3|, particularly the 
relative number of the centrally located rectilinear vor- 
tices and the curved peripheral vortices which have one 
or both ends on the cylindrical side wall. During changes 
of the rotation velocity the curved peripheral vortices 
behave very differently from the central rectilinear vor- 
tices: They have no annihilation energy barrier [3] during 
a deceleration of 51, while during acceleration they may 
become unstable with respect to loop formation and re- 
connection at the cylindrical side wall which leads to the 
generation of new independently evolving vortices in the 
increased rotational flow Jl. 



For simple flow geometries at low velocities vortex fil- 
ament calculations are today a reasonably efficient and 
reliable way to find the resulting vortex configurations, in 
many cases with less effort than by direct measurement. 
The rotating inclined cylinder is a good example since 
an extensive series of measurements would be needed to 
scan the behavior as a function of tilt angle, rotation ve- 
locity, and mutual friction. Here we are not attempting 
to provide the entire quantitative picture of the possi- 
ble flow states, but look for the main qualitative changes 
from the ideal case. 

Already in the 1980's Mathieu et al.^ considered a 
cavity that was tilted with respect to the rotation axis. 
They used the free-energy principle combined with the 
HBVK continuum approximation. They also implicitly 
assumed that far from the boundaries the vortices are 
aligned along the rotation axis, and considered a rect- 
angular cavity which is infinite in one direction (that is 
perpendicular to the rotation axis). In the rectangular 
geometry considered by Mathieu et al.[b\ the superfluid 
velocity without vortices is always along the direction 
which is considered to to be infinitely long. This forces 
the vortices to lie in the plane perpendicular to this di- 
rection. In reality the cavities are not infinitely long and 
often the longest dimension is oriented along the rota- 
tion axis. Therefore the flow of the superfluid component 
without vortices takes a more complicated form. Unless 
the container has rotational symmetry and its symme- 
try axis is aligned perfectly along the rotation axis, the 
superfluid is not at rest any more but moves due to the 
container boundaries. This results in that the superfluid 
velocity is not simply caused by vortices but contains a 
potential part that must be taken into account. 

In the following wc consider a circular cylinder and find 
its vortex configuration in rotation using the vortex fila- 
ment model. The stable vortex configuration is obtained 
by following the evolution in the vortex dynamics long 
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enough in order to obtain the steady state response. The 
most striking result, which differs from the general view 
and from the results assumed in Ref . , is that far from 
boundaries the vortices may tend to align not along the 
rotation axis but along the symmetry axis of the rotating 
cylinder. This is especially true for a single vortex under 
high rotation and for an array with a large number of 
vortices. These configurations would result if the cylin- 
der with a given number of straight vortices is tilted by 
an angle which is less than a critical value. With large 
enough tilt angles the vortices may become unstable and 
tend to orient more along the rotation axis. Above the 
critical angle there exist turbulent regions in the tilted 
cylinder which act as a source for new vortices. There- 
fore the tilt of the cylinder has a similar effect as rotation 
around more than one axisQ. 

The article is organized as follows. First in Sec. [Til we 
solve the velocity profile for the superfluid component 
inside the cylinder when no vortices are present. In the 
following section we describe briefly the vortex filament 
model and the numerical procedure that is used to deter- 
mine the vortex motions and the steady state inside the 
cylinder. Section [IV] is devoted to the results. 



II. 



MOTION OF THE IDEAL FLUID INSIDE 
TILTED ROTATING CYLINDER 




FIG. 1: Cylinder tilted by an angle 77 with respect to the 
rotation axis fl = QZ. 

The helium superfluids, ^He and ■^He, can be described 
using a two-fluid model where the superfluid part behaves 
like an ideal Eulerian fluid with negligible viscosity and 
the normal component behaves like a classical fluid with 
finite viscosity. We are mostly interested in superfluid 
■^He-B, where the normal component is very viscous and 
can always be considered to be in solid body rotation 
in the rotating cavity. If there are no vortices, the su- 
perfluid part is in potential flow, which must be taken 



into account before one can consider the motion of the 
vortices. In case of an ideal circular cylinder, with ra- 
dius R and length L, aligned perfectly along the rotation 
axis, the superfluid remains stationary in the laboratory 
frame. When the cylinder axis is not parallel to the ro- 
tation axis the superfluid is no more at rest but is put 
into motion due to the boundaries. The same occurs also 
in all other cases without cylindrical symmetry, such as 
the rotating boxQ or a cylinder with ellipsoidal cross 
sectionfsl]. Consider a situation (at some particular time) 
illustrated in Fig. [1] where fl = flZ and the cylinder axis 
is along cos?7Z -I- sinryX. At that particular time the 
fixed Cartesian coordinates (A, Y, Z), denoted by upper- 
case letters, and the Cartesian coordinates (x, y, z), fixed 
to the rotating cylinder and denoted by lowercase letters, 
can be related by a rotation by an angle 77 around the Y- 
axis. The cylindrical coordinates (p, 4>, z) that are also 
fixed to the rotating cylinder, such that p — yjx^ + y'^ 
and (f) = arctan(y/a;), can then be related to the fixed 
laboratory coordinates via 



X ~ p cos (j) cos rj + z sin rj 
Y = p sin (j) 

Z = — p cos (jjsmi] + z COST]. 



(1) 



The boundaries of the cylinder are given by the surfaces 
p — R and z — ±L/2. In cylindrical coordinates f2 and 
r are given by 

r2 = r2 — sin 77 cos 0/5 -|- sin 77 sin (/)^ -|- cos 77Z (2) 
V ^ PP + zz. (3) 

Without vortices the superfluid velocity follows that in 
potential flow, v = V<I>, where the potential $ satisfies 
the Laplace equation. The boundary condition is that 
at the surfaces (with normal n) the velocity component 
normal to the boundary satisfies ^ = n • (Jl x r) which 
implies that one must solve the following system of equa- 
tions: 
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The above equations give the solution in the laboratory 
frame. Conversion to the rotating frame is obtained by 
removing the effect of rotation, so that the velocity is 
given by 



V* - f2 X r. 



(5) 



Equation (j4|) can be further simplified with the use of 
a trial function $(p, 0, z) = Q.[pz -\- g{p,z))sni4>smrj, 
which results in that ^(p, z) has to satisfy the following 
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Helmholtz equation: 

dg 
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Note that the potential g{p, z) does not depend on tilt 
angle rj. The velocity components in the rotating frame, 
obtained from Eq. are then given by 



Vp = S2 sm 0sm7i— 
op 

va, — r2cos0sin?7— — ilpcos?7 
P 

Vz = sin(/)sin7y(2p + — ). 

oz 



(6) 



One notices that the tilt of the cylinder implies a flow par- 
allel to the cylinder axis and also reduces the azimuthal 
counterflow. The axial velocity is largest at the cylin- 
der wall (p = i?, z = 0) and varies as sin </) as one goes 
around the cylinder axis. The solution to the above equa- 
tions is found numerically, except in case of an infinitely 
long cylinder, when g = Q. As an example, the solution 
for the vortex-free flow state is plotted in Fig. [2] for the 
potential z) and in Fig. [3] for the velocity v{x, y, z), 
when L/i? = 4and7? = 30°. 




FIG. 2: Solution for velocity potential g{p, z) when L = 47?. 



III. VORTEX FILAMENT MODEL 

Superfluid vortices are modeled as line defects. This 
model is suitable especially for superfluid ^He where the 
core diameter is of the order of atomic size. Our cal- 
culations are attempting to model primarily superfluid 
■^He-B where two simplifications can be assumed: 1) The 
normal fraction is always in solid-body rotation with the 
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FIG. 3: Velocity in the rotating cylinder when L — AR and 
when the angle between the rotation axis and cylinder axis, 
■q — 30°. Left: Velocity in the a:: = plane, the plane that 
is perpendicular to the plane defined by the cylinder and the 
rotation axes. Right top: Velocity in the z = plane. Right 
bottom: Velocity in the bottom of the cylinder, z = —L/2. 
The coordinates {x,y,z) are fixed to the rotating cylinder. 



container, as noted in Sec. [TTl 2) In addition the core 
diameter qq « 10~^ mm. This is still a small value 
compared to the experimental dimensions of a few mil- 
limeters and the use of the vortex filament model is ap- 
propriate. The larger coresize, however, makes it more 
justifiable to assume smooth walls for the cylinder, since 
so far there has been no clear experimental indication 
that vortices become pinned if the walls of the container 
are well polished[9j. Also down to intermediate tempera- 
tures of T > OATc mutual friction damping is sufficiently 
high so that vortices are rapidly annihilated at zero fiow 
conditions. In such conditions vortex remanence is not 
a major problem and relatively large vortex-free normal 
fiuid-superfluid counterflow can be achieved. In practice 
this means that the rotating container can generally be 
filled with any number of vortices up to the number which 
corresponds to the equilibrium state with solid-body ro- 
tation of the superfluid component (lol| . 

The vortex filament model is well described by Schwarz 
in Refs. [Til. Our numerical scheme is similar. Line 
vortices are described by a sequence of points {s^j^Q. 
The required tangent, curvature radius etc. at some par- 
ticular point s(r) along the line vortex are calculated by 
fitting a circle through three points (the point under con- 
sideration and the nearest neighbours). The main differ- 
ence, compared to calculations of Schwarz, is that we use 
a 4th order Runge-Kutta method to determine the vortex 
motion that is described by the following equation: 

VL = Vs + as' X (v„ - Vs) - a's' x [s' x (v„ - Vg)] . (7) 

Here vl = ds/dt is the velocity of the vortex point at 
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s with unit tangent s'. The above equation includes 
the mutual friction between vortex lines and the nor- 
mal fluid with coefficients a and a' that are taken as a 
function of temperature from the measurements of Se- 
van et al. pisl . [ij] This equation is typically used in the 
laboratory frame where Vn — fl xjr, but is also valid in 
the rotating frame where Vn — 0.[l5| Now in the rotating 
frame the superfluid velocity Vg = v -I- v^^ -I- Vb consists of 
the velocity due to rotation, v, calculated above in Sec. 
Inland given by Eq. the superfluid velocity due to vor- 
tices, "v^, and the boundary induced velocity, Vb, that is 
needed to cancel the velocity through the boundaries due 
to vortices. 

As one may notice the above method for doing the cal- 
culation in the rota ting frame differs from previous sim- 
ulations in Refs. 

[3, 

ll7l. [T8t , which contain a complicated 
term, Siot, that was used to take into account the rota- 
tion. This earlier approach cannot be completely correct. 
A simple way to check this is to consider the zero temper- 
ature limit in an infinite system with no other external 
velocities (except the rotation of the normal component, 
which at r = cannot couple to the vortex motion, sim- 
ply because it is absent). In this case our v reduces to 
V — — ft X r and that takes care of the rotating frame. 
This term does not change the vortex configuration, but 
just rotates the coordinates, like expected. However, the 
term s^ot may have a more dramatic change on the vortex 
configuration since the term depends on the vortex con- 
figuration itself. Therefore, owing to its complex form, 
this term has different and configuration-dependent con- 
sequences. 

The superfluid velocity due to vortices is obtained us- 
ing the Biot-Savart formulation: 

K /■ (si - r) X dsi 

Here the integration is a line integral along the vortices 
and K is the quantum of circulation. When r = s is one of 
the vortex points, the above equation becomes singular 
as Si —f s. This singularity can be avoided, for example, 
by considering a moving vortex ring ^llj . If one denotes 
with Z_ and 1+ as being the lengths of two adjacent line 
elements connected to the point s, then the velocity of 
this point is given by 

/■' (Si - s) X rfsi 
4^7 |si-s|3 ■ 

The integral (the term that is typically denoted as the 
non-local term) is now taken over the rest of the vortex 
filament, not connected to the point s. Vectors s' and s" 
in the local term are the tangent and principal normal at 
s and the derivation is with respect to the arc length, ^, 
along the vortex. 

As already described by Schwarz, the boundary veloc- 
ity, Vb = V<I>b, in general, must be obtained by solving 



the Laplace equation V^$b = with the boundary con- 
dition that n-(v^-l-Vb) = on the boundary, with normal 
n. This method has been used in several previous cal- 
culations, for example in Ref. p^ . by using finite differ- 
ences to discretize the Laplace equation and solving the 
resulting sparse matrix equation. The resulting matrix 
is typically very large and limits the grid size that can 
be used to evaluate the boundary velocity. In these cal- 
culations here we used an approximate method, namely 
image vortices to evaluate the boundary velocity. Above 
and below the cylinder we used simple image vortices to 
cancel the velocity through the bottom and top plates of 
the cylinder. The velocity through the cylindrical wall 
at p = i? is canceled using an inverse image curve where 
the distance of the image vortex from the cylinder axis 
is given by pi,„ = R^/p- Here p is the distance of the 
vortex from the axis. The z and coordinates of the 
image vortex are the same as those for the real vortex 
(additionally one needs to reverse the direction of the 
image vortex). A similar method was used by Zieve and 
Donev[20|. This image vortex method is accurate in case 
of straight vortices parallel to the cylinder axis. When 
comparing the simulations with the full solution of the 
Laplace equation, the agreement is quite good since only 
a minor part of the velocity comes from the image field. 

The vortex filament model cannot describe vortex re- 
connection processes and tell when a reconnection occurs. 
We use commonly accepted criteria for making a recon- 
nection: A reconnection in simulations is done when two 
vortices approach each other closer than some small dis- 
tance. This distance, denoted by ft.rcc, is typically taken 
to be about one half of the smallest space resolution used 
to describe the vortices, /imin- Additionally, we require 
that in every reconnection the vortex length must de- 
crease, which approximately indicates that the total en- 
ergy decreases. This additional requirement will prevent, 
for example, two parallel straight vortices from recon- 
necting. The distance between vortices is calculated, not 
simply by determining the distance between the points, 
but by calculating the smallest distance between differ- 
ent vortex segments. These vortex segments (extending 
from Si to Si-|_i are assumed to be straight. Reconnection 
with a cylinder wall is assumed to occur when a vortex 
is closer than /irec/2 from the wall. For plane bound- 
aries this means that the distance to the image vortex 
is smaller than hrec used for bulk reconnections. When 
a reconnection is done, the point closest to the wall is 
split into two and both these points (one being the first 
point on a vortex and the other one being the last point) 
are moved to the surface. The point separation along 
the vortex is kept between some reasonable limits, such 
that ft-niin < — Si| < /iniax and more points are used 
in regions where the vortex is strongly curved. At the 
boundaries the vortex tangent is fixed to be along the 
surface normal. This requirement ensures that there is no 
flow through the surface. Small vortex loops are removed 
when it becomes impossible to follow their dynamics. In 
our numerical scheme this is implemented by removing 
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any closed loop or vortex terminating on the wall which 
is shorter than 'ih^i„. 



IV. RESULTS 

A. Configuration of a single vortex 

Before analysing how a vortex array behaves in the 
tilted cylinder it is informative to consider a single vor- 
tex. The stability of the vortex with respect to creation of 
new vorticesQ depends on mutual friction a(T) but also 
on the initial configuration (which in many experiments 
is unknown). The lower the temperature, the more easily 
a vortex becomes unstable, leading to the creation of a 
new vortex after reconnection at the boundary. All this 
is consistent with experiments P. [l6j . However, with zero 
tilt angle ij only very few configurations lead to multipli- 
cation of vortices, see for example Refs. [3, HH. There- 
fore, to investigate the propagating vortex front [13] and 
the twisted vortex state [l9| behind it, the simulations 
had to be started with a large number of vortices. At 
low temperatures with small mutual friction our simula- 
tions show that the probability for a vortex to become 
unstable increases as the tilt angle increases from zero. 
In order to observe easy vortex multiplication (for vari- 
ous initial configurations) the tilt angle should be much 
larger than the experimental one, which is only of order 
one degree. 

In this article we are not focusing on the criteria for the 
vortex stability with respect to reconnection on the wall 
and the formation of independently moving new loops. 
These criteria seem to depend on many different features 
and parameters (at least on tilt angle, temperature, rota- 
tion velocity and initial vortex configuration) . We are in- 
terested in the (possibly metastable) steady states in the 
cylinder, with different number of vortices in the cylin- 
der. 

The easiest way to find the steady state configurations 
of a single vortex is to start with the vortex along the 
cylinder axis with zero tilt. Then by increasing the tilt 
angle (which was done in Fig. [4] in steps of 10 degrees) 
one may follow the vortex dynamics by iterating Eq. ([7]) 
until the vortex finds its equilibrium position. To speed 
up the process and to avoid vortex multiplication these 
simulations were done at high enough temperatures, with 
the ^He-B parameter values at T = Q.STc- The numerical 
resolution for these single vortex calculations was such 
that /imin — 0.04 mm and /imax = 0.10 mm. The cylinder 
radius was fixed to R = 3 mm. 

For large rotation velocities there was no indication of 
hysteresis, i.e. the obtained configuration was the same 
even if one started the iteration from a straight vortex 
perpendicular the cylinder axis. The steady state con- 
figurations for high rotation velocity are illustrated in 
Fig. 21 One may notice that for a single vortex the steady 
state depends strongly on the aspect ratio L/R of the 
cylinder. For short cylinders L ^ R the vortex is approx- 
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FIG. 4: Configuration of a single vortex in the tilted rotating 
cylinder at a high rotation velocity of f2 = 500 mrad/s and 
radius R = 3 mm. Tilt angles are 77 = 0, 10, 20, . . . , and 
90 degrees. The vortex is practically in the xz-p\a.ne that 
is defined by the rotation and cylinder axes. Top: L — R. 
Bottom left: L = 47?. Bottom right: L = 20i?, only the upper 
half of the vortex configuration is plotted. 



imately along the rotation axis, except near boundaries, 
where the vortex must be along the surface normal. For 
long cylinders, L 3> -R, the vortex lies mainly along the 
cylinder axis and only deviates from this near the bottom 
and top plates. When 77 90° the vortex finally turns 
and aligns itself along the rotation axis. For long cylin- 
ders the vortex is along the cylinder axis because then the 
azimuthal counterflow is reduced in a longer region, even 
if the reduction at every point is not an optimal one. For 
the short cylinder it is more important how well the coun- 
terflow is reduced locally and therefore a better result is 
obtained by orienting the vortex along the rotation axis. 
One should note that, with respect to the number of vor- 
tices, these single vortex states at high rotation velocities 
are metastable and far from equilibrium. However, they 
are experimentally accessible states in superfluid '^He-B, 
where at high temperatures one may easily create either 
a vortex free state or a state with an arbitrary number 



6 



1.5- 



0.5: 



z/R 



-0.5- 



-1.5- 




z/R 




FIG. 5: (Color online) Steady state vortex configurations for 
a single vortex in a tilted cylinder with R = 3 mm and L = 
12 mm. The vortex is in the plane defined by the rotation 
and cylinder axes (a; z- plane). Left: Low rotation velocity of 
10 mrad/s. The tilt angles for vortices terminating to the top 
and bottom ends of the cylinder are: 0, 10, 20, 30, 40, 50, and 
60 degrees and were obtained by increasing the tilt angle. The 
vortices terminating to the outer cell wall (p — R) correspond 
to tilt angles of 0, 10, . . ., 90 and were obtained by decreasing 
the tilt from 90 degrees. The configurations for the highest 
tilt angles (from 70 to 90) can also be obtained by increasing 
the tilt. Right: Rotation velocity of 25 mrad/s. For vortices 
attached to the top and bottom plates of the cylinder the tilt 
angles are 0, 10, 20, . . ., 60 degrees. For vortices terminating 
to the outer cell wall rj = 30, 40, . . ., 90 degrees. Hysteresis 
in increasing/decreasing inclination is observed for angles 30, 
40, 50, and 60 degrees. 



of vortices ^1^. 

At small rotation velocities one may observe hystere- 
sis between different metastable configurations. This is 
quite natural since the boundary condition that a vortex 
must terminate perpendicular to the solid wall will result 
in a self-induced velocity due to its curvature that is of 
order Vcurv = Kln(8i?/ao)/47ri?, where R is the radius of 
the cylinder. This translates to a characteristic rotation 
velocity of order 10 mrad/s at which hysteresis becomes 
important. This is illustrated in Fig.[5]where we have fol- 
lowed the vortex configuration by increasing the tilt angle 
first from zero to 90 degrees and then decreased it back to 
zero (also in steps of 10 degrees). Hysteresis also results 
from the cell geometry and the boundary requirement at 
the solid wall. Therefore a vortex may get trapped either 
on the top/bottom plate or on the cylindrical side wall, 
p = i?, in a similar fashion as a vortex is trapped on a 
spherical pinning site on the plane boundary 11 1 . 



B. Vortex array in an infinitely long cylinder 

When the rotation and cylinder axes are perfectly 
aligned, the vortices form a lattice of rectilinear lines 
which are parallel to the common axis. Closer to the 
cylindrical outer wall the vortex array deforms to con- 
centric rings. These are separated from the cylinder wall 
by a vortex-free annulus which has a thickness of order 
of the intervortex distance P, [1) If one considers an 
infinitely long cylinder, which naturally is impossible to 
organize experimentally, but which gives a good hint of 
what might occur in a long cylinder far from the ends, 
some analytical results are possible. In this case the 
potential velocity profile due to the tilt is obtained by 
putting g — in Eq. One may easily note that 

an array of rectilinear vortices along the cylinder axis 
can still be a stable solution. However, now the areal 
vortex density riy = 2fl cos tj/k is reduced by the factor 
cos r] due to the reduced azimuthal counterflow in Eq. ^ . 
The axial velocity Vz in Eq. (O does not affect the vor- 
tex configuration as long as it is below the critical value 
of the Ostermeier-Glaberson instability [H, H^l- Consid- 
ering the stability of a single vortex and taking into ac- 
count the maximum axial velocity Wz,max = 2i?n sin 77 and 
the reduced counterflow, one obtains that the vortices 
suffer the instability when 

^z,max ^ 2^ COS 7^, where 
J/ = Kln(l/A:ao)/(47r) w k. This results in that the insta- 
bility occurs when 



V cos rj 



sm rj 



(10) 



with the wave number k = cot rj/R. In this derivation we 
assumed that the vortex array fills the whole cylinder and 
we ignored the vortex free region near the walls. Alter- 
natively, by considering a vortex cluster with N vortices, 
with radius Rc = \J kN /2ttQ, cos rj^ placed in the center of 
the cylinder, one obtains that the vortices at the largest 
axial flow (p = R^ and sincj) = ±1) suffer an instability 
when 



iV > iVc = Stt- cot^ w 27r cot^ ry, 



(11) 



which depends only on the tilt angle of the cylinder. 
This is rather well satisfied numerically if one addition- 
ally takes into account that for small N the vortex bun- 
dle is always somewhat smaller than what the analytical 
estimate suggests. Using periodic boundaries along the 
cylinder axis and a period that is much larger than the 
possible wavelength of Kelvin waves, one obtains that 
with a tilt angle rj = 30°, the corresponding Nc — 19 
according to Eq. (jlip . while numerically we find that a 
bundle with 22 vortices is still stable. With 25 vortices 
one may observe growing Kelvin waves, first at locations 
where the axial velocity is largest. 

We have not determined the vortex configurations of an 
infinitely long cylinder in detail, since the computer code 
is not fully compatible with periodic boundaries. How- 
ever, some calculations have been done for finite cylin- 
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FIG. 6: (Color online) Stable vortex array configuration, 
when f2 = 50 mrad/s in a cylinder with L = AR — 12 mm 
and a tilt angle of 15 degrees. Left: Side view along the y- 
axis where the red straight line denotes the rotation axis. Top 
right Top view where the crosses denote the equilibrium lo- 
cation of straight vortices at zero tilt, the configuration used 
as an initial guess for the calculations.. Bottom right: Per- 
spective view of the same configuration. During the itera- 
tion 10 out of the 38 initial vortices annihilate, which indi- 
cates that the counterflow at the cylindrical boundary needs 
to be large enough to prevent the vortices that terminate to 
the outer boundary from shrinking away. The calculations 
were done using mutual friction parameters that correspond 
to T = O.STc for ^He-B. A similar configuration can be ob- 
tained at lower temperatures where it only takes much longer 
time to reach equilibrium, owing to much smaller mutual fric- 
tion. Numerical resolution is characterized by /imin = 0.04 
mm and /imax = 0.25 mm. 



ders, the case more appropriate experimentally. Some of 
these findings are reported below. 



C. Steady states in finite cylinder 

For small rotation velocities and small tilt angles one 
may observe a static steady state. This state is most 
likely not the absolute energy minimum. Frequently the 
calculations for straight vortices have shown that differ- 
ent initial configurations will lead to metastable states 
with very similar energies Some of these steady states 
in the tilted cylinder are shown in Figs. [6] and [T] These 
configurations were obtained from an initial array of 
straight vortices along the cylinder axis as the starting 
state and then iterating the equation of motion with a 
given tilt angle. As in the single vortex case, one may 
notice that in a long cylinder the resulting vortex bundle 
is oriented more along the cylinder axis than along the 
rotation axis. Additionally, the vortices become slightly 




FIG. 7: (Color online) Stable vortex array configuration, 
when f2 = 50 mrad/s in a cylinder with L = lOR = 30 mm 
and a tilt angle of 15 degrees. Left: Side view along the y-axis 
where the red straight line denotes the rotation axis. Center. 
Side view along the x-axis. Top right: Top view. Bottom 
right: Perspective view of the same configuration. During 
the iteration 15 out of the 38 initial vortices annihilate. The 
calculations were done using mutual friction parameters that 
correspond to T = O.STc in "^He-B. The initial configuration 
and the numerical resolution were the same as in Fig. [S] 



twisted, which results from the tendency to reduce the 
axial counterflow. This twist cannot be seen in the case 
of a single vortex, since the vortex lies in the plane where 
the axial flow is absent. 

With large enough tilt and rotation velocity the initial 
array of straight vortices suffers the Glaberson instabil- 
ity, as described above for an infinitely long cylinder. Our 
simulations show that the equation for the critical num- 
ber of vortices, Eq. gives an estimate of the instabil- 
ity also in a cylinder of finite length. If the number of vor- 
tices exceeds this value, then at the location of the largest 
axial flow the vortices go unstable. This is illustrated in 
Fig. [5]where the simulations were started with 30 straight 
vortices at the center of the cylinder. The configura- 
tion at 750s is shown in Fig. [9l together with the time 
development of vortex length and average polarization. 
Pi — / d^- The initial vortex array, which has 

the polarization = 1 and = Py = 0, is broken by the 
Glaberson instability. This has the result that the polar- 
ization along the z-axis decreases at the same time as 
the vortex length increases. The final steady state value 
for p^ = 0.80 is higher than what would result if all the 
vortices would be along the rotation axis, in which case 
pI = cos^ry = 0.75 (and pi = sin^ 77 = 0.25). Since the 
turbulent fluctuations (Kelvin waves along vortices) and 
also the bending of vortices towards the radial direction 
near p = R tend to decrease the polarization along the 
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FIG. 8: (Color online) Development of the dynamic steady 
state, when L = 6R = 18 mm for a rotation velocity of 250 
mrad/s and a tilt angle of 30 degrees. The mutual friction 
parameters correspond to T = OATc in He-B. The initial con- 
figuration is an array of 30 rectilinear vortices that become 
unstable due to the Glaberson instability. Similar break-up 
of the vortex array can also be seen in simulations done at 
T = O.STc- Numerical resolution is characterized by ftmin ~ 
0.10 mm and /imax ~ 0.30 mm. See also Fig. |9l 

2;-axis one may safely argue that even in this dynamical 
steady state there is a tendency of the vortices to orient 
themselves along the symmetry axis of the cylinder. 

A somewhat surprising result from this simulation (and 
some other similar ones which we have performed at 
somewhat different temperature or with different L) is 
that after the instability it is difficult to obtain a static 
steady state, which presumably would be the thermody- 
namic equilibrium configuration. Even after the vortex 
length and polarization of the vortex tangle have settled 
to a steady state value one may still observe creation of 
new vortices in regions of high axial counterflow (which 
arises from the potential flow) owing to the Glaberson in- 
stability and a reconnection with the cylinder boundary. 
These newly generated vortices are then driven towards 
less turbulent regions. At the same time there is contin- 
uous annihilation of small vortex loops near the cylinder 
ends in regions of small counterflow. It is possible that 
our simulations do not extend far enough in time and 
that a static solution follows at much later times. An- 
other possibility is that the calculation gets stuck in a 
local energy minimum and a possible static solution at 
lower energy may not be found with this initial configu- 
ration of vortices. This would also mean that experimen- 
tally it might be possible to observe a similar dynamic 
state. 

At high temperatures the vortices become smoother 
and the turbulent regions, that can be seen in Figs. [5] 
and [5] at T = OAT^ become much more damped, but 
surprisingly some dynamics may still remain. For exam- 




time (s) 



FIG. 9: (Color online) Structure of the tangled bundle of vor- 
tices at 750 s in more detail for simulations with L = 6R = 
18 mm at a rotation velocity of 250 mrad/s and a tilt angle 
of 30 degrees. The mutual friction parameters correspond to 
T = OATc in He-B. Left: View along the y-axis that is perpen- 
dicular to the plane defined by the rotation and cylinder axes. 
Center. View along the s-axis. Top right: Top view along the 
cylinder axis. Bottom right: Time development of the polar- 
ization pi = J (s' • ifd^/ f and the total length JZ — J 
of the vortex tangle. This length £ is scaled by the combined 
total length of rectilinear vortices in an array at zero tilt, in 
the continuum approximation that is {20,/ k)-kFI? L. 



pie, Kelvin waves are still created in the regions of largest 
axial flow. Similarly some vortices shrink in the regions of 
small counterflow near the cylinder ends. A comparison 
between the low and high temperature states is shown in 
Fig.[Tni The low temperature state at T = 0.22Tc was ob- 
tained from an initial unstable configuration with single 
vortex in the form of a quarter ring starting from the bot- 
tom of the cylinder and bending to the outer wall. At low 
temperatures with a tilt angle of 30° this configuration 
is unstable, leading to the creation of new vortices via 
reconnection at the boundary. At high temperatures the 
same initial configuration does not lead to vortex multi- 
plication. Therefore the high temperature state, which 
here corresponds to T r^^Tc (using a = 10 and a' = 0.9), 
was obtained using the state from the T = 0.22Tc calcu- 
lation at t=1500 s as initial configuration and iterating 
the vortex motion with high mutual friction. The steady 
state which was obtained this way was still a dynamical 
one after 300s of evolution, but the creation and annihi- 
lation of vortices was a quite regular process. This can 
be seen, for example, by investigating the vortex length, 
which shows small fluctuations but also clear oscillations 
with a periodicity of order 10 seconds. Fig. fTHl shows also 
that the steady state vortex length increases and the ax- 
ial polarization decreases with decreasing temperatures. 
This is due to fluctuations (Kelvin waves) on the vor- 
tices which are suppressed at high temperatures. The 
difference in steady state vortex length between the low 
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FIG. 10: (Color online) Comparison between the low and high 
temperature vortex configurations for a short cylinder when 
L = 5 mm, R — 3 mm, Q — 0.25 rad/s and i] = 30°. Now 
the vortices are mostly along the rotation axis (noted by the 
red line in the uppermost plots). From top to bottom these 
three different views correspond to a view along the y-, x-, 
and a-axis, respectively. The low temperature configuration 
(taken at t = 1800 s) on the left is for T = 0.22Te in ^He-B 
and was obtained by using a single vortex, a quarter ring, as 
an initial configuration. The high temperature configuration 
on the right, which corresponds to mutual friction parameters 
a = 10 and a' — 0.9, was obtained by using the configuration 
at t = 1500 s from the simulations at T = 0.22Tc as an initial 
guess and iterating for 262 seconds. Bottom left: Time de- 
velopment of the total vortex length, scaled by {2Q/ KjnR^ L. 
The solid black line corresponds to T = 0.22rc and the solid 
red line, which starts at t= 1500 s, is the high temperature 
limit. The dashed blue line is a zero temperature calcula- 
tion where we have set a = a' = and the only dissipation 
is due to numerics. The initial configuration for the T 
calculation was the same as for T = 0.22Tc. Bottom right: 
Time development of the polarization of the vortex tangle. 
Color-coding is the same as for the vortex length. Large os- 
cillations at small times result from the motion of a single (or 
few) vortex around the cylinder. 



and high temperature states is of order 20%. In reality 
the difference is larger since simulations cannot describe 
structures that are smaller than the numerical resolu- 
tion, which in these runs is determined by hy^in ~ 0.10 
mm and /imax = 0.30 mm. A small overshoot in vortex 
length, that can be seen in simulations at low temper- 
atures before the length saturates, indicates that a sur- 
plus of vortices is created before the final steady state is 
reached. Further analysis of the present simulations and 
more calculations are required to understand the com- 
plex dynamic phenomena in a tilted cylinder at large tilt 
angles with a large number of vortices. For example, 
the verification of Eqs. (|10p and pT|) is far from com- 
plete, especially when fie (or equally when A^c) is large. 
Unfortunately the simulations take a lot of time and to 
determine whether the steady state is static or dynamical 
one is a tedious task. 



V. CONCLUSIONS 

We have investigated vortex motions in a rotating 
cylinder where the rotation axis is not parallel to the 
symmetry axis of the cylinder. We have determined the 
steady state solutions and obtained that for small tilt 
angles the stable configuration consists of vortices that 
are along the cylinder axis, rather than the rotation axis. 
Therefore the assumptions behind the previously calcu- 
lated result might be inappropriate although there the 
calculations concerned a different fiow geometry. For 
large tilt angles (but still well below 90°) the steady 
state solution, especially at low temperatures, may be 
a dynamical one, which would make it possible to study 
a polarized turbulence in a well controlled environment, 
especially in superfluid '^He-B where strong vortex pin- 
ning can be ignored. Since the vortex configuration is 
different from the expected one, it might be appropriate 
to reconsider the idea of a tilt-dependent mutual fric- 
tion parameter, B(r\), as used in Ref. Q to interpret 
second-sound measurements in superfluid ^He. However, 
one should note that in He-II pinning is important and 
might alter the results. 



Acknowledgments 

I would like to thank M. Krusius for his valuable com- 
ments and improvements to the article. This work is 
supported by the Academy of Finland (grant 114887). 



10 



[1] D. Stauffer and A.L. Fetter, Phys. Rev. 168, 156 (1968). 
[2] J.A. Northby and R.J. Donnelly, Phys. Rev. Lett. 25, 
214 (1970). 

[3] V.M. Ruutu, J.J. Ruohio, M.Krusius, B. Plagais, and 

E.B. Sonin, Pliysica B 255, 27 (1998). 
[4] A. P. Finnc, V.B. Eltsov, G. Eska, R. Hanninen, J. Kopu, 

M. Krusius, E. V. Tliuncberg, and M. Tsubota, Phys. 

Rev. Lett. 96, 085301 (2006). 
[5] P. Mathieu, B. Plagais, and Y. Simon, Phys. Rev. B 29, 

2489 (1984). 

[6] M. Kobayashi and M. Tsubota, Phys. Rev. A 76, 045603 

(2007); J. Low Temp. Phys. 150, 587 (2008). 
[7] L.M. Milne-Thomson, Theoretical Hydrodynamics 5th ed. 

(Macmillan, New York, 1968). 
[8] L.D. Landau and E.M. Lifshitz, Fluid Mechanics (Perg- 

amon Press, Oxford, 1963). 
[9] V.B. Eltsov, R. de Graaf, R. Hanninen, M. Krusius, R.E. 
Solntsev, V.S. L'vov, A.I. Golov, and P.M. Walmsley, in 
Progress in Low Temperature Physics: Quantum Turbu- 
lence, Vol. 16, 45-146 (2009); 
[10] V.M.H. Ruutu, U. Parts, J.H. Koivuniemi, N.B. Kopnin, 

and M. Krusius, J. Low Temp. Phys. 107. 93 (1997). 
[11] K.W. Schwarz, Phys. Rev. B 31, 5782 (1985). 
[12] K.W. Schwarz, Phys. Rev. B 38, 2398 (1988). 
[13] T.D.C. Bevan, A.J. Manninen, J.B. Cook, H. AUes, J.R. 
Hook, and H.E. Hall, J. Low Temp. Phys. 109, 423 
(1997). 



[14] T.D.C. Bcvan, A.J. Manninen, J.B. Cook, A.J. Arm- 
strong, J.R. Hook, and H.E. Hall, Phys. Rev. Lett. 74, 
750 (1995). 

[15] E.B. Sonin, Rev. Mod. Phys. 59, 87 (1987). 

[16] A. P. Finnc, T. Araki, R. Blaauwgccrs, V.B. Eltsov, N.B. 

Kopnin, M. Krusius, L. Skrbck, M. Tsubota, and G.E. 

Volovik^ Nature 424, 1022 (2003). 
[17] M. Tsubota, C.F. Barenghi, T. Araki, and A. Mitani, 

Phys. Rev. B 69, 134515 (2004). 
[18] M. Tsubota, T. Araki, and C.F. Barenghi, Phys. Rev. 

Lett. 90, 205301 (2003). 
[19] V.B. Eltsov, A. P. Finnc, R. Hanninen, J. Kopu, M. Kru- 
sius, M. Tsubota, and E. V. Thuneberg, Phys. Rev. Lett. 

96, 215302 (2006). 
[20] R.J. Zieve and L.A.K. Donev, J. Low Temp. Phys 121, 

199 (2000). 

[21] R. Hanninen, A. Mitani, and M. Tsubota, J. Low Temp. 
Phys. 138, 589 (2005). 

[22] V.B. Eltsov, A.I. Golov, R. do Graaf, R. Hanninen, M. 
Krusius, V.S. L'vov, and R.E. Solntsev, Phys. Rev. Lett. 
99, 265301 (2007). 

[23] W.I. Glaberson, W.W. Johnson, and R.M. Ostermeier, 
Phys. Rev. Lett. 33, 1197 (1974); R.M. Ostermeier and 
W.I. Glaberson, J. Low Temp. Phys. 21, 191 (1975). 

[24] R.J. Donnelly, Quantized Vortices in Helium 11 (Cam- 
bridge University Press, Cambridge, UK, 1991). 



